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Abstract 

A new gregarious braconid parasitoid wasp of Euphorinae, Meteorus stellatus Fujie, Shimizu & Maeto 
sp. nov., is described from the Ryukyu Islands in Japan, based on an integrative taxonomic framework. 
The phylogenetic position of the new species within the Meteorini was analyzed based on DNA frag- 
ments of the mitochondrial cytochrome c oxidase 1 (CO1) and the nuclear 28S rDNA genes. The new 
species was recovered as a member of the versicolor complex of the versicolor + rubens subclade within the 
pulchricornis clade. The new species is a gregarious parasitoid of two Macroglossum species (Lepidoptera: 
Sphingidae) and constructs single or several unique star-shaped cocoon masses separately suspended by 


very long threads. The evolution of gregariousness and spherical cocoon masses is discussed. 
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Introduction 


The pupae of parasitoid wasps cannot actively escape various risks, such as predation, 
parasitism, pathogenesis, and environmental stresses. Therefore, cocoons and mum- 
mies play important roles in protecting soft and exarate pupae from such risks (Gauld 
and Bolton 1988; Shaw and Huddleston 1991; Maeto 2018). 

Members of Braconidae, one of the most diverse hymenopteran families, form 
various types of cocoons and mummies to adapt to various natural enemies and en- 
vironmental threats. Many gregarious braconids produce communal cocoon masses, 
while discrete cocoons for each individual are also constructed (Shaw and Huddles- 
ton 1991). Final-instar braconid larvae emerge from host organisms and then spin a 
cocoon using silk from the labial glands (e.g., Shaw and Huddleston 1991; Quicke 
2015). In some gregarious species, wasp larvae not only spin individual cocoons but 
also cooperate to construct massive communal cocoons (Zitani and Shaw 2002). 

The cosmopolitan braconid genus Meteorus Haliday consists of more than 300 
valid species (Yu et al. 2016). As the genus Ze/e Curtis is apparently nested within the 
Meteorus species tree, Meteorus is a paraphyletic group and its rearrangement into sev- 
eral monophyletic genera is pending (Maeto 1990b; Stigenberg and Ronquist 2011). 
Maeto (1990b) has divided the Meteorus species into seven morphologically and bio- 
logically defined species groups, while monophyly has been supported partially (Sti- 
genberg and Ronquist 2011). 

Meteorus species are solitary or gregarious koinobiont endoparasitoids of Lepidop- 
tera or Coleoptera larvae (Huddleston 1980; Stigenberg and Ronquist 2011; Maeto 
2018). Their cocoons are either suspended (pendulous) or not, and the cocoon suspen- 
sion is one of diagnostic characters to identify the pulchricornis group (Maeto 1989a). 
The gregarious species construct either independent (e.g., /. acerbiavorus Belokobyl- 
skij, Stigenberg & Vikberg, M. heliophilus Fischer, and M. rubens Nees) or communal 
cocoons (e.g., MZ. congregatus Muesebeck, M. komensis Wilkinson, M. kurokoi Maeto, 
and M. townsendi Muesebeck) (Maeto 1989b; Zitani and Shaw 2002; Zitani 2003; 
Stigenberg et al. 2011). The latest molecular phylogeny has suggested that the gregari- 
ousness is a derived character in the pulchricornis + rubens group complex (Stibenberg 
and Ronquist 2011). So far, six types of cocoon architectures have been observed in 
the gregarious species of Meteorus, as is shown in Table 1. However, the evolution of 
gregariousness and cocoon architectures within the genus Meteorus has been poorly 
studied, and more bionomical and phylogenetic information is thus needed. 

Although the cocoon structure of Meteorus is quite mysterious, larval behavior as- 
sociated with cocoon formation has received little attention, with only a few reports ex- 
amining it (Askari et al. 1977; Zitani and Shaw 2002; Maeto 2018). Recently, unique 
star-shaped cocoon masses of an undescribed gregarious species of Meteorus have been 
observed in subtropical Japan (Mitamura 2013). Therefore, this study aims to describe 
that gregarious species of Meteorus based on integrative morphological and molecular 
evidence and observe its cocoon-mass formation behavior. The phylogenetic position of 
the new species and evolution of gregariousness and cocoon masses are also discussed. 
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Table |. The types of cocoon masses in gregarious Meteorus. 


Types Characteristics 


A loosely clumped within a host pupal chamber 


B individually suspended from host plant by a thread 

C sparsely arranged and suspended by a common cable 

D loosely clumped and suspended by a common cable 

E congregated and directly attached to host plant without a cable 
F communal and suspended by a common cable 


Species 
Palearctic species 
M. acerbiavorus 
M. heliophilus 
M. rubens 
Neotropical species 
M. oviedoi 
M. papiliovorus 
Palearctic species 
M. kurokoi 
Neotropical species 
M. restionis 
Neotropical species 
M. cecavorum 
M. juliae 
Neotropical species 
M. congregatus 
Afrotropical species 
M. komensis 
Neotropical species 
M. townsendi 
undescribed species 
Oriental species 
M. stellatus sp. nov. 
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Present study 


Materials and methods 


Study fields 


The field collection of host moth larvae to observe the cocoon formation behavior of 
emerged larvae of wasps was conducted at Okinawa Municipal Museum, Okinawa 
City, Okinawa-honté, Okinawa Prefecture, Japan. Some materials were also collected 
within Okinawa-honté (Okinawa Prefecture) and Amami-déshima (Kagoshima Pre- 
fecture), Japan. All materials were from the middle part of the Ryukyu Islands, the 
subtropical Oriental region in Japan. 


Morphological observation and terms 


Morphological observation was conducted with a stereoscopic microscope (SMZ800N, 
Nikon, Tokyo, Japan). Specimens and cocoons were photographed using a Digital 
Microscope (VHX-1000, Keyence, Osaka, Japan) with a 10-130x lens. Multi-focus 
photographs were stacked in the software associated with the Keyence System. Mullti- 
focus photographs of cocoon masses were taken using a single lens reflex camera («71, 
Sony, Tokyo, Japan) with a micro-lens (A FE 50 mm F2.8 Macro SEL50M28, Sony). 
The RAW format photographs were developed using Adobe Lightroom CC v.2.2.1 
(Adobe Systems Inc., San Jose, CA, USA), and stacked using Zerene Stacker v.1.04 
(Zerene Systems LLC., Richland, WA, USA). The holotype of MZ. komensis, deposited 


in the Natural History Museum, London, UK was also examined by the second author 
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using a stereoscopic microscope (SMZ1500, Nikon). Multi-focus photographs were 
taken using an «7II camera with micro-lenses (LAOWA 25 mm F2.8 2.5—5x ULTRA 
MACRO, Anhui Changgeng Optics Technology Co., Ltd, Hefei, China). The cap- 
tured RAW format photographs were developed and stacked as per the aforementioned 
photo technique used for the cocoon masses. The figures were edited in Microsoft 
PowerPoint 2019. 

The description style mostly follows that of Stigenberg and Ronquist (2011). The 
morphological terms and measurements follow those of Richards (1977) and van 
Achterberg (1988). The following abbreviations are used: OOL = ocelli-ocular line, 
OD = ocelli diameter of a posterior ocellus, and POL = posterior ocellar line. 

The abbreviations for repositories are listed below: 


CNC Canadian National Collection of Insects, Ottawa, Canada; 

ELKU Entomological Laboratory, Faculty of Agriculture, Kyushu University, 
Fukuoka, Japan; 

EMUS Utah State University Insect Collection (= American Entomological Institute: 
AET), Department of Biology, Utah State University, Logan, Utah, USA; 

EUM Ehime University Museum, Matsuyama, Japan; 

KPMNH_ Kanagawa Prefectural Museum of Natural History, Odawara, Japan; 

MNHA = Museum of Nature and Human Activities, Sanda, Japan; 

NARO Institute for Agro-Environmental Sciences, NARO (= NIAES: National 
Institute for Agro-Environmental Sciences), Tsukuba, Japan; 

NHMUK_ Natural History Museum, London, United Kingdom (formerly BMNH); 

NSMT National Museum of Nature and Science, Tsukuba, Japan; 

OMM Okinawa Municipal Museum, Okinawa, Japan; 

OMNH = Osaka Museum of Natural History, Osaka, Japan; 

RUM Ryukyu University Museum, Okinawa, Japan; 

SEHU the Laboratory of Systematic Entomology, Faculty of Agriculture, Hok- 
kaido University, Sapporo, Japan; 

TARI Taiwan Agricultural Research Institute Council of Agriculture, Executive 
Yuan, Taichung, Taiwan; 

ZISP Zoological Institute, Russian Academy of Sciences, St Petersburg, Russia. 


Secondary sex ratio 


To investigate the secondary sex ratio of wasps, the number of males and females of all 
enclosed wasps that emerged from each host larva was counted for 11 host larvae (Sup- 
pl. material 1: Table $1). Cocoon masses with hyper-parasitoids were excluded, and the 
number of dead cocoons was not counted. The effect of the number of total wasps on 
the proportion of males was analyzed using a generalized linear model (GLM). A logit 
link function and binomial error distribution were employed. The mean proportion 
of males was estimated from the average number of total wasps. ‘The analysis was per- 


formed in IBM SPSS Statistics for Windows, v. 25.0 (IBM Corp., Armonk, NY, USA). 
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Observation of cocoon formation behavior 


The cocoon formation behavior of wasp larvae was observed at a laboratory of OMM, 
in June 2019 by the third author. It was recorded with video cameras (Sony Handy- 
cam, HDR-CX470 and HDR-XR150, Sony). The suspended larvae were blown with 
air currents created by breathing, as there was no wind, which would enhance the 
merging of each individually suspended larva in the laboratory as it would in natural 
conditions. A single silk thread spun by an individual larva is called a "thread", and 
intertwined threads are called a "cable" as in Barrantes et al. (2011). A short movie 
showing the cocoon formation behavior is available on YouTube (https://www.you- 


tube.com/watch?v=AuHarLHolPM). 


Molecular analysis 


Gene selection 


To delimit a species, fragments of a mitochondrial protein encoding gene, cytochrome 
c oxidase 1 (CO1), were selected, because its evolutionary rate is more or less rapid and 
it is one of the most common genes used for population to species level phylogenetic 
analysis (this is well-known as the DNA barcoding gene). To infer the phylogenetic 
relationships among species of Meteorini (Meteorus and Zele), CO1 and a nuclear 
noncoding gene, 28S rRNA (28S), were selected. 28S is a gene that has evolved more 
slowly than CO1 and is usually used for species-groups or higher-level phylogeny; 
therefore, the combined CO1 and 28S analysis can provide a higher resolution of spe- 
cies phylogeny. 


Taxon sampling and outgroups 


A total of 44 species of Meteorus including five morphospecies were sampled as in- 
groups. Five species of Ze/e were also sampled as ingroups because Zele is deeply nested 
within the Meteorus tree (Stigenberg and Ronquist 2011). Three species from differ- 
ent tribes of euphorine genera were sampled as outgroups (Syrrhizus Forster, Syntretus 
Forster, and Peristenus Forster). A total of 177 sequences of CO1 and 172 of 28S 
were compiled from GenBank. Sequences obtained from databases sometimes contain 
unreliable information (e.g., Klimov et al. 2019; Shimizu et al. 2020), and an evalu- 
ation of such sequences is always strongly recommended to ensure that the analysis is 
accurate. In the present study, the sequences used by Stigenberg and Ronquist (2011) 
and several additional sequences were considered as reliable sources. The complete 
information of the sampled taxa and sequences is available in the Suppl. material 2: 


Table S2. 


DNA extraction, amplification, and sequencing 
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Table 2. Primer information for PCR. 


Target Primer name Sequence (5’ to 3”) References 
Col CO1 Ico hym CAA ATC ATA AAG ATA TTG G Schulmeister (2003) 
CO1 hco extB* CCT ATT GAW ARA ACA TAR TGA AAA TG Schulmeister (2003) 
28S 28SD1F ACC CGC TGA ATT TAA GCA TAT Harry et al. (1997) 
28SD5R* CCC ACA GCG CCA GTT CTG CTT ACC Schulmeister (2003) 


The newly collected samples from Okinawa were stored in 99.9% ethanol for DNA 
extraction. DNA was extracted from a right mid or/and hind leg. The protocols followed 
from PCR to sequencing were according to the work of Shimizu et al. (2020), except for 
the primers in Table 2, for which the PCR conditions were as follows: CO1: initial dena- 
turation (2’00”) at 95 °C, 35 cycles of denaturation (0’30”) at 95 °C, annealing (0°30”) at 
48 °C, an extension (1’00”) at 72 °C, and a final extension (10’00”) at 72 °C; 28S: 1°00” 
at 95 °C, 40 cycles of 0°30” at 95 °C, 0’30” at 48 °C, 1°30” at 72 °C, and 5’00” at 72 °C. 


Species delimitation 


Partial fragments of CO1 were used for species delineation. A total of 189 sequences were 
used for analysis (Suppl. material 2: Table $2). To obtain an accurate multiple sequence 
alignment (MSA), MSA was conducted using MEGA v.10.0.5 (Kumar et al. 2018) based 
on amino acids. First, the codon positions of all nucleotide sequences were adjusted, the 
nucleotide sequences were translated to amino acids, the amino acid sequences were 
aligned by CLUSTAL W (Thompson et al. 1994) implemented in MEGA with default 
settings, the amino acid alignment was checked by eye, and finally the aligned amino acid 
sequences were retranslated to nucleotides. ‘The final dataset was 657 bp without indels. 

To delimit the species, both distance- and topology-based methods were employed 
as below. Using both methods, three types of datasets were analyzed: (1) Meteorus + 
Zele + outgroups, (2) Meteorus + Zele, and (3) Meteorus. Prior to the analysis, identical 
haplotypes were removed from the datasets on the web server of ALTER (Glez-Pefa et 
al. 2010) (available at: http://sing.ei-uvigo.es/ALTER/). 


Distance-based method (ABGD) 


The barcoding gap based analysis, Automatic Barcode Gap Discovery (ABGD) (Puil- 
landre et al. 2012), was ran on the graphic web version of ABGD (available at: http:// 
wwwabi.snv.jussieu.fr/public/abgd/abgdweb.html), with the following parameters: 
Pmin = 0.001, Pmax = 0.1, Steps = 10, X (relative gap width) = 1.0, model = Kimura 
(K80) (TS / TV = 2.0), and Nb bins (for distance distribution) = 20. 


Topology-based method (GMYC) 
The General Mixed Yule Coalescent (GMYC) analysis was employed. GMYC analysis 


requires an ultrametric tree (UTtee) as an input. To construct the UTree, the model and 
parameters was selected on a web server of the smart model selection (SMS) (Lefort et al. 
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2017) (available at: http://www.atgc-montpellier.fr/phyml-sms/): the GTR+G+I model 
was selected as the best fit model under the Bayesian information criterion (BIC). The 
UTtee was generated using BEAST v.2.6.3 (Bouckaert et al. 2019), with a random start- 
ing tree, the uncorrelated lognormal relaxed clock model, and the coalescent tree prior. 
A Bayesian Markov chain Monte Carlo (MCMC) was run for 40,000,000 generations, 
with trees sampled every 5,000 generations, and a burnin of anterior 25%. The con- 
vergence of run was assessed using Tracer v.1.6 (Rambaut and Drummond 2007): run 
reached a stationary distribution and all effective sample sizes (ESS) were greater than 
200. A majority-rule consensus ultrametric tree was finally generated using TreeAnnota- 
tor v.2.6.3 (Bouckaert et al. 2019). GMYC analysis was run using the GMYC function 
of the R-package Splits (Fujisawa and Barraclough 2013) (available at: http://r-forge.r- 
project.org/projects/splits/) using R v.3.6.3 (R Core Development Team 2020). 


Phylogenetic analysis 


The phylogeny of Meteorini species was inferred with both the Bayesian Inference 
(BI) and maximum likelihood (ML) approaches using a concatenated CO1 and 28S 
fragments. 


MSA 


Although the MSA of CO1 was already performed in the species delimitation, MSA 
for 28S was conducted in the MAFFT online service (Katoh et al. 2019), using the 
Q-INS-i algorithm, which is the structural alignment method for RNA (Katoh and 
Toh 2008; Katoh and Standley 2013). Ambiguously aligned regions were automati- 
cally removed from the dataset using trimAl v.1.2 (Capella-Gutierrez et al. 2009), with 
default parameters. The final datasets were 657 (CO1) and 560 (28S) bp in length: the 
concatenated CO1 and 28S dataset was 1,217 bp. 


Terminal species selection 


In order to exclude the taxon sampling bias, a single sequence for each species was se- 
lected based on the conservative results of the species delimitation analysis by ABGD: 
sequences of 61 Meteorus species, six Zele species, and three outgroup species were 


finally selected (Table 3). 


Model selection 


Each codon position within the CO1 fragment was treated as a different data block, 
but not for noncoding 28S. The best-fit substitution model was determined using Par- 
titionFinder v.2.1.1 (Lanfear et al. 2017) with the greedy search algorithm under the 
corrected Akaike information criterion (AICc): the selected model was the GTR+I+T 
model for the first and second codon position of CO1 and 28S, and GTR+I for the 
third codon position of CO1. 
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Table 3. Nomenclature systems for Meteorini species. The following abbreviations are used: pulchri. = 


pulchricornis. 
Species 


Meteorus ictericus A, B 
M. ruficeps 

M. aff. ruficeps 
Meteorus sp. 

M. artocercus 

M. cinctellus 

M. colon 

M. stenomastax 
M. tenellus A-C 
M. pendulus A-C 
M. abscissus 

M. limbatus 

M. pulchricornis 
Meteorus sp. 

M. acerbiavorus 
M. rubens A-C 
M. aff. versicolor A, B 
M. arizonensis 

M. obsoletus 

M. tarius 

M. stellatus sp. nov. 
M. versicolor 

M. micropterus 

M. abdominator A, B 
M. affinis A-D 

M. cespitator A, B 
M. cis A, B 

M. consimilis 

M. densipilosus 

M. eklundi 

M. filator A, B 

M. gigas 

M. hirsutipes 

M. jaculator 

M. kyushuensis 

M. longicaudis 

M. obfuscatus 

M. aff. obfuscatus 
M. oculatus 

M. sibyllae 
Meteorus sp. 
Meteorus sp. 

M. sulcatus 

M. tabidus 

M. vexator 

Zela albiditarsus 
Z. caligatus 

Z. chlorophthalmus 
Z. deceptor 

Z. niveitarsis 


Zela sp. 


Analysis 


Present study Stigenberg et al. (2011) 
clade / subclade / complex Clade 

ictericus | — | — I 

ictericus | — | — I 

ictericus | — | — I 


ictericus | — | — = 


pulchri. / colon | — HA 
pulchri. / colon | — HA 
pulchri. / colon | — HA 
pulchri. / colon | — HA 
pulchri. / colon | — HA 
pulchri. / pendulus | — HB 
pulchri. / pulchri. /— IB 
pulchri. / pulchri. /— IIB 
pulchri. / pulchri. /— IIB 
pulchri. / pulchri. /— - 
pulchri. / rubens-versicolor | rubens HB 
pulchri. / rubens-versicolor | rubens IIB 


pulchri. / rubens-versicolor | versicolor - 
pulchri. / rubens-versicolor | versicolor - 
pulchri. / rubens-versicolor | versicolor HB 
pulchri. / rubens-versicolor | versicolor - 
pulchri. / rubens-versicolor / versicolor - 


pulchri. / rubens-versicolor | versicolor HB 
micropterus | — | — IIC 
Unresolved Ill 
Unresolved Ill 
Unresolved Ill 
Unresolved Ill 
Unresolved Ill 
Unresolved Ill 
Unresolved Ill 
Unresolved Ill 
Unresolved = 
Unresolved Ill 
Unresolved Ill 
Unresolved Ill 
Unresolved Ill 
Unresolved Ill 
Unresolved Ill 
Unresolved Ill 
Unresolved Ill 
Unresolved 7 
Unresolved ote 
Unresolved Ill 
Unresolved Ill 
Unresolved Ill 
Zele IV 
Zele IV 
Zele IV 
Zele IV 
Zele a 
Zele - 


Maeto (1990) 
group / subgroup 
ictericus | — 
ictericus | — 
ictericus | — 


pulchri. / colon 
pulchri. / colon 


pulchri. / gyrator 
pulchri. / gyrator 
pulchri. / pulchri. 


rubens | — 
pulchri. / versicolor 
pulchri. / versicolor 


pulchri. / versicolor 
micropterus | — 


hirsutipes | — 
hirsutipes | — 


hirsutipes | — 


hirsutipes | — 


The BI analyses were conducted using MrBayes v.3.2.2 (Ronquist et al. 2012). A Bayes- 
ian MCMC analysis was ran with the following settings: four independent runs, 20 
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chains each, heating 0.05, random starting trees, and trees sampled every 1,000" gen- 
eration for 10,000,000 generations. The convergence of the MCMC runs was checked 
by the average standard deviation of split frequencies (ASDSF) in MrBayes (i.e., AS- 
DSF < 0.01) (Ronquist and Huelsenbeck 2003) and chain stationarity in Tracer v.1.6 
(Rambaut and Drummond 2007). Then, we discarded the anterior 25% of the genera- 
tions as burn-in, obtained estimates for the harmonic means of the likelihood scores 
from the remaining 75% of the generations using the sump command, and conducted 
a final check of the convergence of the runs by the value of a potential scale reduction 
factor (PSRF); if the runs were convergent enough, PSRF was less than 5% divergent 
from 1.0. Finally, a consensus tree with the Bayesian inference posterior probabilities 
was obtained using the sumt command in MrBayes. The ML analysis was conducted 
in IQ-TREE v.2.1.2 (Minh et al. 2020) with a Shimodaira-Hasegawa-like approximate 
likelihood ratio test (SH-aLRT) (Guindon et al. 2010) and ultrafast likelihood boot- 
strap replicates (UFBoot2) (Minh et al. 2013; Hoang et al. 2018) for 10,000 replicates. 
The trees were edited in Fig Tree v.1.4.3 (Rambaut 2006-2016), Adobe Illustrator CC 
v.23.0.2, and Photoshop CC v.20.0.4 (Adobe Systems Inc., San Jose, CA, USA). 


Results 


Species identification 


Although the molecular species delimitation was conducted using (1) whole data- 
sets (i.e., Meteorus plus Zele plus outgroups), (2) Meteorus plus Zele datasets, and (3) 
Meteorus datasets, the results were congruent among all datasets in both the ABGD and 
GMYC methods. The number of recognized species was higher in the GMYC than in 
the ABGD method (Fig. 1). However, all results indicated that Meteorus stellatus sp. 
nov. was a single species. 

Baded on morphological data, M. stellatus sp. nov. ran to the versicolor subgroup of 
the pulchricornis group based on Maeto’s (1989a, b) criteria, but it could not be identi- 
fied as any described species of Meteorus (see also Differential diagnosis). 


Taxonomic account 


Meteorus stellatus Fujie, Shimizu & Maeto, sp. nov. 
http://zoobank.org/D8785F79-E874-4854-95D7-5C0A9289 14CA 
[Japanese name: Hoshigata-haraboso-komayubachil] 


Figs 2, 3 


Etymology. The specific name is a masculine Latin word, “stellatus”, meaning “starry”, 
which is derived from the unique shape of the cocoon masses. 

Type series. 41992 40¢'¢ (all from Japan). Holotype 9 (OMNH): “Japan: Ok- 
inawa Is., Okinawa City, Goya / Kitanakagusuku Vil., Shimabukuro, Okinawa Ko- 
domono-kuni, 3. VI. 2019 cocoon masses), K. TONE et al. leg.” “9. VI. 2019 emerged”. 
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ABGD GMYC 
MN NNN KJ591525 Syntretus zuijleni 


MN MEN K 591511 Streblocera carinata 
EEE EE MG 926869 Peristenus sp. 
MN NNN KJ591465 Meteorus sp. 

MN MEN HQ263962 Meteorus aff. ruficeps 


— peerage | Meteorus ruficeps 
HQ263958 . 


MN MEN HQ263956 Meteorus ictericus A 


HQ263957 

HQ264016 ure 
HOQ263928 Meteorus ictericus B 
HQ264013 


HQ263930 
HQ263960 
MEE MEE HQ263973 Zele caligatus 


HQ263938 ] 
( i? Hazese6e | 22/2 chlorophthalmus 


MN MEN HQ263965 Zele albiditarsus 
MEE MEN KJ591540 Zele sp. 

NN, ENN KJ591541 Zele niveitarsis 
HQ263937 
HQ263964. 


t HQ263926 
HQ264028 
Haz64032 Meteorus consimilis 
HQ264033 


MN MEN K.591543 Meteorus sp. JS10-00532 
CJ Ha264025 
[as 
[cas 
HQ263976 


C3 HO263939 | 7 
= az64026 Meteorus affinis B 


HQ263977 
HQ263978 


HQ264051 
HQ263922 | Meteorus affinis D 
HQ263975 


MN NEN HQ263913 Meteorus tabidus 
(EE MEE HOQ263921 Meteorus abdominator A 
MN EN HQ.263954 Meteorus abdominator B 
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Figure |. Species delimitation of Meteorus plus Zele plus closely related outgroups based on ABGD and 


GMYC methods, shown using a Bayesian consensus ultrametric tree generated using BEAST. Although 


the species delimitation was conducted using (1) Meteorus plus Zele plus outgroups, (2) Meteorus plus Zele, 


and (3) Meteorus sequences or topologies, all results were congruent; therefore, all results are shown as a 


summarized unit for each method. 
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Paratypes: 12233 (OMM), same as holotype; 22 9254 (OMM), Okinawa- 
kodomonokuni, Goya, Okinawa City /Shimabukuro, Kitanakagusuku Vil, Okinawa- 
honté Is., collected as cocoon masses on 13.V.2019 and emerged on 19.V.2019, Koichi 
Tone et. al leg; 2292¢¢ (OMNH), Nago City, Okinawa-honté Is., collected as 
cocoon masses on Morinda umbellata and emerged on 23—24.V.2011, Masashi Sugi- 
moto leg.; 49 24¢0¢ (OMNB), Nago City, Okinawa-honté Is., collected as host larva 
of Macroglossum passalus passalus, and adult wasps emerged on XII.2010 (22 9233) 
and 1.2011 (29924), Masashi Sugimoto leg.; 229243 (RUM), Hentona, 
Kunigami Vil., Okinawa-honté Is., collected as cocoon masses on 14.IV.2011 and 
emerged on 19.IV.2011, Koichi Sugino leg.; 29 02¢¢ (ZISP), Hentona, Kunigami 
Vil., Okinawa-honté Is., collected as cocoon masses on 30.X.2010 and emerged on 
7.X1.2010, Kozue Miyagi leg.; 29 O2¢¢ (NARO), Aha, Kunigami Vil., Okinawa- 
honté Is., collected as cocoon masses on 23.X.2010 and emerged on 25.X.2010, 
Yasuji Kakazu leg.; 29 0233 (OMNH), Okinawa-kodomonokuni, Goya, Okinawa 
City, Okinawa-honté Is., collected as host larva of Macroglossum pyrrhosticta feeding 
Paederia scandens on 3.V1.2019, cocoon masses formed on 9.V1.2019, and emerged 
on 17.VI.2019, Yu Erh Chen leg.; 29 92¢.¢ (CNC), Ryukyu University, Nishihara 
Town, Okinawa-honté Is., collected as cocoon masses on 22.V1.2005 and emerged on 
25-29.V1.2005, Kazuo Minato leg.; 202 2203.3 (ELKU, EMUS, EUM, KPMNH, 
MNHA, NHMUK, NSMT, OMM, SEHU and TARI with 22 9244 each), Yoshi- 
hara, Chatan Town, Okinawa-honto Is., collected as cocoon masses on 4. VI.2020 and 
emerged on 9-11.VI.2020, Tamami Gushiken leg.; 12 (OMNH), Chuo-rindo, Ama- 
mi City, Amami-éshima Is., 5. VII.2013, Shunpei Fujie leg. 

Non-types. 3239 92284'¢ adults; 29 cocoon masses (see Suppl. material 1: Table S1). 

Distribution. Japan (Ryukyus: Okinawa-honté Island and Amami-éshima Island). 

Differential diagnosis. Meteorus stellatus sp. nov. is most similar to VM. komensis 
(Fig. 4) but can be distinguished from the latter by the following combination of char- 
acters: the comparatively larger posterior ocelli (in M. stellatus sp. nov., OOL/OD = 
1.2—1.6, while 1.7 in MZ. komensis (Fig. 4C)), the medially longitudinally strigose first 
metasomal tergite (Figs 2J, 3F) (fairly striate in /. komensis, as is shown in Fig. 4E, F). 
Meteorus stellatus sp. nov. is also very similar to M. kurokoi but can be distinguished 
from M. kurokoi by the face width (1.5—1.8x its height in M. stedlatus sp. nov. (Figs 2E, 
3B), while 1.2-1.3x in M. kurokoi) and the sculpture of frons (frons with a median 
longitudinal or a pair of carinae in MW. stellatus sp. nov. (Fig. 2G), but smooth without 
any carinae in VM. kurokoi). 

In the key to species of Meteorus from the West Palaearctic region (Stigenberg and 
Ronquist 2011) and China (Chen et al. 2004), MZ. stellatus sp. nov. would run to M. 
versicolor (Wesmael), but can be distinguished from M. versicolor by the face width (1.5— 
1.8x its height in M. stellatus sp. nov. (Figs 2E, 3B), while 1.0-1.2x in M. versicolor), the 
sculpture of frons (with a median longitudinal or a pair of carinae in M. stellatus sp. nov. 
(Fig. 2G), but smooth in M. versicolor), and the shape of temple (Figs 2C, 3C) (roundly 
narrowed in /. stellatus sp. nov. (Figs 2C, 3C), but directly narrowed in M. versicolor). 
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J . ‘ , + = ™ 2 j . L wy pene 
Figure 2. Meteorus stellatus sp. nov., 2 holotype (exceptionally L is a paratype) A habitus B head, frontal 


view C head, dorsal view D mesopleuron and scutellum, dorsal view E mesosoma, lateral view F head, 
dorso-lateral view G forewing H basal antennal segments I apical antennal segments J propodeum and 


T1, dorsal view K T2 and following tergites, dorsal view L T1, ventral view. 


The results of a GenBank BLAST search showed that the CO1 sequences of 
M. stellatus sp. nov. were closest to those of M. arizonensis Muesebeck and M. tarius 
Huddleston. However, M. stellatus sp. nov. can be distinguished from M. arizonensis 
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by its smaller body (the body lengths of M. stellatus sp. nov. and M. arizonensis are 
2.9-3.9 mm and 4.6—5.5 mm, respectively), the longer malar space (the malar space 
length 1.0-1.4x the basal mandibular width in /. stellatus sp. nov. whereas 0.6—0.7x 
in M. arizonensis), the shorter ovipositor sheaths (the ovipositor sheath length 1.1—1.2x 
length of the first tergite in M. stellatus sp. nov. and 1.6—-1.9x in M. arizonensis); the 
species can be distinguished from M. tarius by its broader face (the face width 1.5—1.7x 
its height in . stellatus sp. nov. whereas approximately 1.0x in M. tarius) and the posi- 
tion of the forewing vein m-cu (slightly antefurcal to interstitial in MW. stellatus sp. nov., 
but far antefurcal in VM. tarius). 

Description. Female (holotype; Fig. 2). Body length 3.6 mm. 

Head (Fig. 2B, C, E, H, I). Width of head 1.7x median height. Length of eye 1.7x 
length of temple in dorsal view. Temple roundly narrowed posteriorly. Eyes large and 
moderately convergent ventrally. Face with width 1.6x height; distinctly and densely 
transversely striate with fine granulation. Clypeus as wide as face, distinctly separated 
from face, and punctate-rugose. OOL / OD = 1.2. POL / OD = 1.4. Frons widely 
smooth, anteriorly with a pair of obscure carinae. Vertex and temple almost smooth. 
Length of malar space 1.1x basal mandibular width. Antennae with 26 segments; 4" 
segment 3.1x longer than wide; and penultimate one 1.9x longer than wide. 

Mesosoma (Fig. 2D, E, J). Mesosoma length 1.4x height. Mesoscutum entirely 
covered with short and dense pale setae; median lobe weakly punctulate in anterior 
0.6 and mostly rugose reticulate in posterior 0.4; lateral lobes weakly punctulate. No- 
tauli shallow, wide, complete, coarsely rugose-reticulate. Prescutellar depression deep, 
almost straight, with often five rather fine carinae. Scutellum smooth and distinctly 
convex. Mesopleuron mostly punctulate, and rugose reticulate anterodorsally. Precoxal 
sulcus shallow and widely rugose-reticulate. Propodeum entirely coarsely rugose re- 
ticulate without median longitudinal carina. 

Wings (Fig. 2G). Fore wing with 3.1 mm in length, length of pterostigma 2.9x 
maximum width, 3-SR/r = 1.2, m-cu distinctly (left) to slightly (right) postfurcal, cu-a 
far postfurcal, 1-CU1 / cu-a = 0.8. Hind wing with 1M / cu-a = 0.8, 1M / Ir-m = 0.6. 

Legs. Tarsal claws with a distinct submedial lobe. Hind leg: outer surface of coxa 
punctate; femur 4.7x longer than wide, and distinctly and densely punctate. 

Metasoma (Fig. 2J, K, L). 1* tergite 1.6x longer than apical width; dorsopes 
absent; mostly smooth anteriorly, longitudinally strigose with some rugosity medi- 
ally, densely striate in posterior 0.3; ventral borders jointed from the base of segment 
to about middle point. Remaining terga smooth. Ovipositor slightly down-curved; 
length of ovipositor sheath 0.7x C+SC+R, 0.3x fore wing, and 1.2x 1* tergite. 

Color (Fig. 2A). Brownish yellow, except for following parts infuscate: stem- 
maticum, apical segments of antennae, dorsal part of propleuron, side of scutellum, 
metanotum, propodeum, apex of hind femur and tibia, posterior half of 1* tergite, 2"4 
tergite (often except for anteromedially), 3 tergite and ovipositor sheath; and palpi 
pale yellow. Wing membrane hyaline; pterostigma light brown, faintly paler basally. 

Variation. Body length 2.9-3.9 mm. Width of head 1.6-1.8x median height. 
Length of eye 1.5—1.7x length of temple in dorsal view. Face with width 1.5—1.7x 
height. OOL / OD = 1.2-1.6. POL / OD = 1.3-1.7. Frons with a longitudinal carina 
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dorsal view D basal antennal segments E apical antennal segments F propodeum and T1, dorsal view 


G T2 and following tergites, dorsal view. 


or a pair of obscure carinae. Length of malar space 1.0—1.4x basal mandibular width. 
Antennae with 26-31 segments; 4° segment 2.9-3.6x longer than wide; and penul- 
timate one 1.7—2.0x longer than wide. Mesosoma length 1.4—1.5x height. Fore wing 
length 2.7—3.5 mm with length of pterostigma 2.8-3.3x maximum width, 3-SR/ r = 
0.8—1.4, m-cu distinctly postfurcal to interstitial, 1-CU1 / cu-a = 0.6-1.1. Hind wing 
with 1M / cu-a = 0.6-1.0, 1M / Ir-m = 0.5-0.7. Hind femur 4.6—4.9x longer than 
wide. 1 metasomal tergite 1.5—1.8x longer than apical width; longitudinally strigose 
with often some rugosity medially; length of ovipositor sheath 0.6—0.8x C+SC+R and 
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Figure 4. Meteorus komensis Wilkinson, 8 holotype A habitus B head, frontal view C head, dorsal view 
D mesoscutum, dorso-lateral view E T1, dorsal view F T2 and following tergites, dorsal view G forewing 


H mesosoma, lateral view. 


1.1-1.2x 1* tergite. 2°¢ tergite brownish yellow to infuscate anteromedially. Pterostig- 
ma unicolored or faintly paler basally. 

Males (Fig. 3). Similar to females, except for length of eye 1.6—1.9x length of tem- 
ple in dorsal view; width of face 1.6-1.8x height; OOL/ OD = 1.1-1.4; POL/ OD = 
1.41.8; length of malar space 1.1—1.7x basal mandibular width; antennae with 27-31 
segments; penultimate segments of antenna 2.2—2.8x longer than wide; hind femur 
4.6-5.1x longer than wide; 1* metasomal tergite 1.6—-1.9x long than apical width; 
propleuron not darkened dorsally. 

Hosts and host plants. Two species of Sphingidae (Lepidoptera) were iden- 
tified as hosts of M. stellatus sp. nov.: Macroglossum passalus (Drury) feeding on 
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Figure 5. Seasonal changes in the adult emergence of Meteorus stellatus sp. nov. (bars indicating the 


number of cocoon masses) and in the monthly average temperature in Naha City, Okinawa-honté Island 


(a solid line). 


Daphniphyllum glaucescens Blume (Daphniphyllaceae) and M. pyrrhosticta Butler feed- 
ing on Paederia foetida Linnaeus [= P scandens (Lour.) Merr.] (Rubiaceae). All wasp 
larvae of M. stellatus sp. nov. emerged from mature larvae of the host sphingids. 
Hyper-parasitoids. Some hymenopteran hyper-parasitoids emerged from the cocoon 
masses after the emergence of /. stellatus sp. nov. adults. The following three species were 
identified as morphospecies at the generic level: Tetrastichus sp. (Eulophidae), Eurytoma 
sp. (Eurytomidae), and Aphanogmus sp. (Ceraphronidae) (Suppl. material 2: Table S2). 
Habitats. Despite the multiple field collection sessions at primary forest areas in 
the Okinawa-honté and Amami-éshima Islands, only one specimen of M. stellatus sp. 
nov. was sampled from a secondary evergreen forest in the latter island. Most other 
specimens of MM. stellatus sp. nov. were collected from a campus of the University of the 
Ryukyus, urban parks, and back yards in Okinawa-honté Island, by finding suspended 
cocoon masses or rearing host larvae. As the host sphingids and their host plants are 
abundant in or around the edges of sparse forests, M. stellatus sp. nov. likely prefers 
rather open forests. 
Phenology. The emergence of adult wasps occurred from April to June and from 
October to January, but not during the hottest season from July to September (Fig. 5). 
Secondary sex ratio. The proportion of males (secondary sex ratio) ranged from 
0.20 to 0.64, showing a gradual increase with the total number of wasps per host 
larva (Fig. 6). The positive effect of the total number of wasps on the sexual ratio was 
significant (B = 0.010, Wald Chi-Square = 18.129, df = 1, p < 0.001). The estimated 
mean of the proportion of males was 0.36 (95% confidence interval = 0.31—0.42, 
significantly less than 0.50) in the average number of total wasps (50.5), indicating an 
overall female-biased sex ratio. 
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Figure 6. Sex ratios of emerged adults of Meteorus stellatus sp. nov. in relation to the number of individu- 


als per host larva. 


Cocoon mass formation. The third author observed a case of larval emergence and 
subsequent cocoon formation in the laboratory. At approximately 1:30 p.m. on June 9, 
2019, approximately 100 larvae of M. stellatus sp. nov. emerged from the abdomen of a 
matured larva of M. pyrrhosticta on the vine of P. foetida by chewing holes (Fig. 7A). The 
larvae commenced hanging down from the host plant substance with their own suspenso- 
ry threads as soon as they emerged (Fig. 7B). The larvae began to cluster by actively rotat- 
ing, twisting, swaying, and horizontally stretching (Fig. 7C). When the larvae found the 
threads of other individuals, they actively went down, intertwined with said threads, and 
eventually merged together. Once they formed a large mass, the mass did not descend any 
more (Fig. 7D). In rare cases, several larvae moved from the upper mass to the lower mass 
as the cable of the former intertwined with that of the latter, owing to the blowing wind. 
No further larval transfer was observed after approximately 70 min of the emergence of 
the larvae. Initially, the shape of the larval masses was irregular (Fig. 7E), but gradually the 
larvae adopted a spherical shape (Fig. 7E, G). The larvae twisted their upper bodies and 
spun the thread at the posterior of their body, namely inside the cocoon mass. ‘The silk 
walling action lasted approximately 40 min, along with the spinning of their own indi- 
vidual cocoons (Fig. 7H, I). Finally, three cocoon masses were completed approximately 2 
h after larval emergence (Fig. 7J). A video of the entire process of cocoon mass formation 
is available at the following address: https://www.youtube.com/watch?v=AuHarLHolPM. 

The host sphingid died on the following day after wondering. The color of the 
cocoons gradually darkened over a few days. After 8 days, 68 females and 23 males 
of M. stellatus sp. nov. emerged from these three cocoon masses. The wasps emerged 
simultaneously, cutting the tip of each cocoon. 
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Figure 7. Cocoon forming behavior of Meteorus stellatus sp. nov. A emerging from a host larva (start 
time) B hanging down from the host plant substance (2 min) C intertwining with threads: arrows show 
larvae looking for other threads (39 min) D almost merging into three masses (57 min) E-J forming 
spherical cocoon masses (E 30 min F 65 min G 69 min H 84 min I 105 min J 139 min). 


Characteristics of the cocoon masses. The cocoon masses of M. stellatus sp. nov. 
(Fig. 8) were light brown to brown, 7—14 mm in width, 9-23 mm in length, and regu- 
larly spherical to ovoid with minimally 12 (Fig. 8D) to maximally over 100 cocoons 
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(Fig. 8C). Exceptionally, approximately 200 cocoons formed a collapsed large mass in 
an artificial breeding case (Fig. 8B). Each cocoon mass was suspended by a single thick 
cable. The cable was 12-100 cm in length. Although most larvae constructed such 
cocoon masses, sometimes a few single larvae formed their own cocoons on the cable 
(Fig. 8F). The cable consisted of individual threads, which were tightly intertwined, 
like a rope (Fig. 8G). The anterior third to half of individual cocoons was exposed 
outward and fairly distributed on the spherical or ovoid surface. The posterior half of 
individual cocoons was invisible under the dense silk wall. Adults emerged by opening 
an anterior outside cocoon cap, which was circular shaped and tapering (Fig. 8F). Co- 
coons with such a regular cap are typical of the pulchricornis clade of Meteorus (Askari 
et al. 1977; Maeto 1989a, b, 1990a). 

Phylogeny of Meteorini and affinity of M. stellatus sp. nov. The Meteorini phy- 
logeny is illustrated in Fig. 9. The Meteorini classification was also revised based on 
our phylogeny and those of Maeto (1990b) and Stigenberg et al. (2011) (Table 3). 
Although our topology was poorly resolved at the species level, it was mostly congru- 
ent with that of Stigenberg and Ronquist (2011). Meteorini was recovered as a mono- 
phyletic group. Meteorini species were divided into five clades (Fig. 9; Table 3). Zele 
was recovered as a robustly supported monophyletic clade and nested within Meteorus 
species. Monophyly was robustly supported for the ictericus and pulchricornis clades 
but not for the unresolved clade. The pulchricornis clade was divided into four internal 
subclades (the colon, pendulus, pulchricornis, and rubens subclades). 

Meteorus stellatus sp. nov. was recovered as an ingroup of the versicolor complex of 
the rubens subclade within the pulchricornis clade and sister to M. tarius. 


Discussion 


Our observation of M. stellatus sp. nov. shows that gregarious cocoon masses were con- 
structed by the highly elaborated cooperation of larvae. The larvae never merged immedi- 
ately after emergence from their host, but initially just descended. ‘This seems to reinforce 
the idea that the suspended cocoon makes the pupating wasp inaccessible to some poten- 
tial enemies (Shaw and Huddleston 1991; Quicke et al. 2006; Zitani 2003; Shirai and 
Maeto 2009; Maeto 2018). It seems that predators, like ants, seldom encounter suspended 
larvae because the threads of larvae are attached to the plant subtract only by a small area. 

The cable of gregarious Meteorus is thought to be very resistant to breaking and highly 
tolerant to environmental stress (Barrantes et al. 2011). A cable of M. stellatus sp. nov. 
consists of a lot of individual threads and seems to be very strong, pretty much like that 
of M. restionis Shaw & Jones (Barrantes et al. 2011). Interestingly, cocoon masses of M. 
townsendi are suspended by a fairly long cable, the length of which is approximately 3.0 m 
(Zitani and Shaw 2002; Zitani 2003), while the longest cable of M. stellatus sp. nov. is ap- 
proximately 1.0 m. According to our observation, the long individual threads and the ad- 
equate wind during the hanging period seem to make it easier for the larvae of M. stellatus 
sp. nov. to merge with each other. ‘The larvae continue to spin thread until they are able to 
merge into a cocoon mass, and in most cases, they never go down unnecessarily after that. 
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Figure 8. Cocoon masses of Meteorus stellatus sp. nov. A habitus, medium-sized B habitus, exceptionally 


large-sized and somewhat collapsed in an artificial condition C a medium-sized cocoon mass D, E small- 


sized cocoon masses F independent cocoons near a cocoon mass G a part of suspending thread, consisting 


of individual cable. 
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Figure 9. Maximum likelihood tree of Meteorini generated using IQ- TREE (BIPP, Bayesian inference 
posterior probabilities; SH-aLRT, a Shimodaira-Hasegawa-like approximate likelihood ratio test; UF- 
Boot2, ultrafast likelihood bootstrap replicates). 


The star-shaped cocoon masses of /. stellatus sp. nov. can reduce the risk of hyper- 
parasitism, because the exposed area of each individual cocoon is apparently smaller 
than the solitary cocoon or non-star-shaped cocoon masses, as suggested by the spheri- 
cal cocoons of M. komensis (Zitani 2003). The outer cocoons of non-circular cocoon 
masses of gregarious Cotesia glomerata (Linnaeus) (Braconidae, Microgastrinae) are 
actually more easily parasitized than the inner ones (Tagawa and Fukushima 1993; 
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Tanaka and Ohsaki 2006). Therefore, the evolution of gregariousness and spherical 
cocoon masses seems reasonable. 

The sex ratio has been studied in gregarious species of Meteorus, while a simi- 
lar pattern of female-biased sex ratio has been shown in Macrostomion sumatranum 
(Enderlein) (Braconidae, Rogadinae), which is also a gregarious parasitoid of matured 
sphingid larvae (Maeto and Arakaki 2005). Both in the gregarious parasitoid category, 
the number of wasps emerged from each host varies widely (8-122 in M. stellatus sp. 
nov. and 26-160 in Ma. sumatranum) and the proportion of males increases with 
it. The female-biased sex ratio could be a result of the local mate competition, as ex- 
pected in inbreeding gregarious parasitoids (Hamilton 1967; Godfray 1994; Smart 
and Mayhew 2009), in which the increase of the male proportion may be caused by 
the oviposition of multiple females on a single host larva (Werren 1983). This predic- 
tion will be tested by the examination of mating systems, oviposition behavior, and 
primary sex ratio. 

The pulchricornis and rubens species-groups belong to the monophyletic line- 
age of pulchricornis clade (Fig. 10), in which both solitary and gregarious species 
are included and cocoons are usually suspended by a spun thread (Maeto 1990b; 
Stigenberg and Ronquist 2011). It is thus likely that communal cocoon masses of 
type F have evolved through individually or sparsely suspended gregarious cocoons 
of type B and C, and subsequently loosely clumped and suspended gregarious co- 
coons of type D, from suspended solitary cocoons within the versicolor complex. 
However, the evolutionary pathways are not yet clarified because only a few gre- 
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Figure 10. Phylogenetic relationships among the species of pulchricornis clade and their lifestyle. 
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garious species are placed in the present phylogram. Further and comprehensive 
analyses including more gregarious species are necessary to confirm and expand this 
evolutionary scenario. 
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